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Abstract 

We present the possibility of calculating the quark flavor changing neutral current decays B® — > tl and 
B® — > ££ for a large variety of supersymmetric models. For this purpose, the complete one-loop calculation 
has been implemented in a generic form in the Mathematica package SARAH. This information is used by 
SARAH to generate Fortran source code for SPheno for a numerical evaluation of these processes in a given 
model. We comment also on the possibility to use this setup for non-supersymmetric models. 



1. Introduction 

With the recent discovery of a bosonic resonance [2 [2] showing all the characteristics of the SM Higgs 
boson a long search might soon come to a successful end. In contrast there are no hints for a signal of 
supersymmetric (SUSY) particles or particles predicted by any other extension of the standard model (SM) 
[3H7] . Therefore, large areas of the parameter space of the simplest SUSY models are excluded. The allowed 
mass spectra as well as the best fit mass values to the data are pushed to higher and higher values jSj . This 
has lead to an increasing interest in the study of SUSY models which provide new features. For instance, 
models with broken i?-parity [9] [10] or compressed spectra [11] might be able to hide much better at the 
LHC, while for other models high mass spectra are a much more natural feature than this is the case in the 
minimal-supersymmetric standard model (MSSM) |12) . 

However, bounds on the masses and couplings of beyond the SM (BSM) models follow not only from direct 
searches at colliders. New particles also have an impact on SM processes via virtual quantum corrections, 
leading in many instances to sizable deviations from the SM expectations. This holds in particular for the 
anomalous magnetic moment of the muon |13] and processes which are highly suppressed in the SM. The 
latter are mainly lepton flavor violating (LFV) or decays involving quark flavor changing neutral currents 
(qFCNC). While the prediction of LFV decays in the SM is many orders of magnitude below the experimental 
sensitivity [2], qFCNC is experimentally well established. For instance, the observed rate of b — > S7 is in 
good agreement with the SM expectation and this observable has put for several years strong constraints 
on qFCNCs beyond the SM [15]. 

The experiments at the LHC have reached now a sensitivity to test also the SM prediction for BR(i?° — > 
Hp) as well as BR(B° up) [TB] 

BR(B° S -> up) SM = (3.23 ± 0.27) • 1(T 9 , (1) 
BR(B% -> up) SM = (1-07 ± 0.10) • 10~ 10 . (2) 
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Recently, LHCb reported the first evidence for B® —> fi[i. The observed rate [T7] 

BR(B° S -> tf) = (3.2±i:|) x 10-° (3) 

fits nicely to the SM prediction. For BR(£>° — > /i/i) the current upper bound: 9.4 • 10 -10 is already of the 
same order as the SM expectation. 

This leads to new constraints for BSM models and each model has to be confronted with these mea- 
surements. So far, there exist several public tools which can calculate BR(£?° d — > ££) as well as other 
observables in the context of the MSSM or partially also for the next-to-minimal supersymmetric standard 
model (NMSSM) Q2]: superiso Q!], SUSY_Flavor gU], NMSSM-Tools [21], MicrOmegas [22] or SPheno [25] . 
However, for more complicated SUSY models none of the available tools provides the possibility to calculate 
these decays easily. This gap is now closed by the interplay of the Mathematica package SARAH [21] and the 
spectrum generator SPheno. SARAH already has many SUSY models incorporated but allows also an easy 
and efficient implementation of new models. For all of these models SARAH can generate new modules for 
SPheno for a comprehensive numerical evaluation. This functionality is extended, as described in this paper, 
by a full one-loop calculation of B® d — >• 11. 

The rest of the paper is organized as follows: in sec. [2] we recall briefly the analytical calculation for 
BR(J3° d ->■ li). In sec. [3] we discuss the implementation of this calculation in SARAH and SPheno before we 
conclude in sec. [4] The appendix contains more information about the calculation and generic results for 
the amplitudes. 



2. Calculation of BR(B s ° d -)• It) 

In the SM this decay was first calculated in ref [25], in the analogous context of kaons. The higher order 
corrections were first presented in [26] ; see also [27) . In the context of supersymmetry this was considered 
in [25] . See also the interesting correlation between BR(_B" — > /j,p,) and (g — 2) M [29 . 

We present briefly the main steps of the calculation of BR(£?° —> Ikii) with q = s,d. We follow closely 
the notation of ref. |30j . The effective Hamiltonian can be parametrized by 

(CsxyOsxy + CvxyOvxy + CtxOtx) , (4) 

X.Y=L,R 

with the Wilson coefficients Csxy,Cvxy,Ctx corresponding to the scalar, vector and tensor operators 
Osxy = {qjPxq i )(eiP Y ek) , Ovxy = {q 3 r P xqi){ill»PYtk) , Otx = (qj^PxtuXfoM ■ (5) 

Pl and Pr are the projection operators on left respectively right handed states. The expectation value of 
the axial vector matrix element is defined as 

(0\h»l 5 q\B° g (p)) = ijffBO. (6) 

Here, we introduced the meson decay constants /go which can be obtained from lattice QCD simulations 
[31] . The current values for B® and B d are given by [32] 

/ B o = (227 ± 8) MeV , f B o = (190 ± 8) MeV . (7) 

Since the momentum p of the meson is the only four-vector available, the matrix element in eq. ^ can only 
depend on p M . The incoming momenta of the b antiquark and the s (or d) quark are p± , P2 respectively, where 
p = Pi + P2- Contracting eq. Q with p^ and using the equations of motion bf^ = —brrib and f^q = m q q 
leads to an expression for the pseudoscalar current 



The vector and scalar currents vanish 



(0 \brfq\ B n q (p)) = (0 \bq\ Blip)) = . (9) 



From eqs. ^ and ^ we obtain 



(0 \h»PL,R.q\ Blip)) = T-V/j, , (0 \bP L/R q\ Bl{p)) = ±l ™ B f*" 9 ■ (10) 

In general, the matrix element M. is a function of the form factors Fs, Fp, Fv,F A of the scalar, pseudoscalar, 
vector and axial-vector current and can be expressed by 

(4tt) 2 X = F s ££ + Fpl^t + Fypyfcfl + F A p^ 1 b £ . (11) 

Note that there is no way of building an antisymmetric 2-tensor out of just one vector p^. The matrix 
element of the tensor operator Otx must therefore vanish. The form factors can be expressed by linear 
combinations of the Wilson coefficients of eq. Q 

i M B»fB° 

Fs = ~. 7 ipSLL + Cslr — Csrr — Csrl) , (12) 

4 m b + m q 

i M| g /flO 

Fp — 7 (— C S ll + C S lr — Csrr + C S rl) , (13) 

4 nib + m q 

1 

Fv = —-^fB°(CvLL + CvLR—CvRR — CvRL), (14) 

1 

Fa = —^7Ib°{—Cvll + Cvlr—Cvrr + Cvrl)- (15) 

The main task is to calculate the different Wilson coefficients for a given model. These Wilson coefficients 
receive at the one-loop level contributions from various wave, penguin and box diagrams, see Figures |B.4f 
B.10 in Appendix B Furthermore, in some models these decays could also happen already at tree-level 



[33]. The amplitudes for all possible, generic diagrams which can contribute to the Wilson coefficients have 
been calculated with FeynArts/FormCalc [34] and the results are listed in Appendix B This calculation has 
been performed in the DR scheme and 't Hooft gauge. How these results are used together with SARAH and 
SPheno to get numerical results will be discussed in the next section. 
After the calculation of the form factors, the squared amplitude is 

(4tt) 4 \M\ 2 - 2 |F S | 2 (M|o - {mi + m fc ) 2 ) + 2 \F P \ 2 (m|o - (m t - m fc ) 2 ) (16) 
+ 2 \F V \ 2 (M 2 B o{m k - m e ) 2 - (to 2 , - m 2 ) 2 ) 
+ 2 \F A \ 2 (M 2 B o{m k + to,) 2 - (to 2 - to 2 ) 2 ) 
+ 43t{F s F^){m e - m k ) (m%o + (m fe + to £ ) 2 ) 
+ m{F P FX)(m e + m k ) (m 2 B o - (m k - m e ) 2 ^j . 

Here, mi and m k are the lepton masses. In the case k = £, this expression simplifies to 

2Ml , . 
\M\ 2 = j^^(\F s \ 2 + \Fp + 2 m iF A \ 2 ). (17) 

Note, the result is independent of the form factor Fy in this limit. In the SM the leading one-loop contri- 
butions proceed via the exchange of virtual gauge bosons. They are thus helicity suppressed. Furthermore, 
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m B o = 5.36677 GeV 
m„o = 5.27958 GeV 



f B o = 227(8) MeV 
f B l = 190(8) MeV 



t b o = 1.466(31) ps 
t b I = 1.519(7) ps 



Table 1: Hadronic input values by default used for the numerical evaluation of B® d — > 11 in SPheno. 



since these are flavor changing neutral currents, they are GIM suppressed. The diagrams involving virtual 
Higgs bosons are suppressed due to small Yukawa couplings. In BSM scenarios these suppressions can be 
absent. 

The branching ratio is then given by 



BR(fl?->44) = 



t b° \M\' 



16tt M b o ^ 

with t b o as the decay constant of the mesons. 



1 - 



m k + mi 
M B o 



\ 



1 - 



m k - mi 



(18) 



3. Automatized calculation of B° . — > ££ 

3.1. Implementation in SARAH and SPheno 

SARAH is the first 'spectrum-generator-generator' on the market which means that it can generate Fortran 
source for SPheno to obtain a full-fledged spectrum generator for models beyond the MSSM. The main 
features of a SPheno module written by SARAH are a precise mass spectrum calculation based on two-loop 
renormalization group equations (RGEs) and a full one-loop calculation of the mass spectrum. Two-loop 
results known for the MSSM can be included. Furthermore, also the decays of SUSY and Higgs particles are 
calculated as well as observables like £i — > £jj, li — > 3£j, b — > sj, 5p, (g — 2), or electric dipole moments. For 
more information about the interplay between SARAH and SPheno we refer the interested reader to Ref . [35 . 



Here we extend the list of observables by BR(£^ 



and BR(£?2 -> U). For this purpose, the 



generic tree-level and 1-loop amplitudes calculated with FeynArts/FormCalc given in Appendix B have 
been implemented in SARAH. When SARAH generates the output for SPheno it checks for all possible field 
combinations which can populate the generic diagrams in the given model. This information is then used to 
generate Fortran code for a numerical evaluation of all of these diagrams. The amplitudes are then combined 



to the Wilson coefficients which again are used to calculate the form factors eqs. (12)-(15). The branching 
ratio is finally calculated by using eq. (18). Note, the known 2-loop QCD corrections of Ref. [55] are not 
included in this calculation. 

The Wilson coefficients for B® d 11 are calculated at a scale Q = 160 GeV by all modules generated by 
SARAH, as this is done by default by SPheno in the MSSM. Hence, as input parameters for the calculation the 
running masses and couplings at this scale are taken. Furthermore, the default input values for the hadronic 
parameters given in Table [T] are used. These are hard-coded at the moment and have to be changed in the 
file LowEnergy_ [MODEL]. f 90, but a support of the Flavor Les Houches Accord (FLHA) [5S] to change these 
values via the input files should follow in the future. Other input parameters which enter the calculation 
like the CKM matrix, the masses of the SM particles or Gf can be changed in the Les Houches input file 
[37] , While SPheno includes the chiral resummation for the MSSM, this is not taken into account in the 
routines generated by SARAH because of its large model dependence. 

3.2. Generating and running the source code 

We describe briefly the main steps necessary to generate and run the SPheno code for a given model: 
after starting Mathematica and loading SARAH it is just necessary to evaluate the demanded model and call 
the function to generate the SPheno source code. For instance, to get a SPheno module for the B-L-SSM, 
use 
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«[$SARAH-Directory] /SARAH. m; 
Start ["BLSSM"] ; 
MakeSPheno [] ; 

MakeSPheno [] calculates first all necessary information (i.e. vertices, mass matrices, tadpole equations, 
RGEs, self-energies) and then exports this information to Fortran code and writes all necessary auxiliary 
functions needed to compile the code together with SPheno. The entire output is saved in the directory 

[$SARAH-Directory]/Output/BLSSM/EWSB/SPheno/ 

The content of this directory has to be copied into a new subdirectory of SPheno called BLSSM and afterwards 
the code can be compiled: 

cp [$SARAH-Directory] /Output/BLSSM/EWSB/SPheno/* [$SPheno-Directory] /BLSSM/ 
cd [$SPheno-Directory] 
make Model=BLSSM 

This creates a new binary SPhenoBLSSM in the directory bin of SPheno. To run the spectrum calculation 
a file called LesHouches. in. BLSSM containing all input parameters in the Les Houches format has to be 
provided. SARAH writes also a template for such a file which has been copied with the other files to /BLSSM. 
This example can be evaluated via 

. /bin/SPhenoBLSSM BLSSM/LesHouches . in . BLSSM 

and the output is written to SPheno . spc . BLSSM. This file contains all information like the masses, mass ma- 
trices, decay widths and branching ratios, and observables: in the entries 43-48 of the block SPhenoLowEnergy 
are the results given for BR(B° -> e+e") (43), BR{B° -> ^+ >~) (44), BR(B° -> fi^e*) (45), BR(B° d -> 
e+e-) (46), BR(B° d -> fi + (J,~) (47), BR(B° -> r+r") (48). 

The same steps can be repeated for any other model implemented in SARAH, or the SUSY-Toolbox scripts 
[3 5) can be used for an automatic implementation of new models in SPheno as well as in other tools based 
on the SARAH output. 

3.3. Checks 

We have performed several cross checks of the code generated by SARAH: the first, trivial check has 
been that in the case of heavy SUSY spectra the SM limits are recovered. Secondly, as mentioned in the 
introduction there are several codes which calculate these decays for the MSSM or NMSSM. A detailed 
comparison of all of these codes is beyond the scope of the presentation here and will be presented elsewhere 
|38j . However, a few comments are in order: the code generated by SARAH as well as most other codes 
usually show the same behavior. There are differences in the numerical values calculated by the programs 
because of different values for the SM inputs. For instance, there is an especially strong dependence on the 
value of the electroweak mixing angle and, of course, of the hadronic parameters used in the calculation 
|16j . In addition, these processes are implemented with different accuracy in different tools: the treatment 
of NLO QCD corrections [35], chiral resummation [JDJ, or SUSY box diagrams is not the same. Therefore, 
we depict in Fig. [T]a comparison between SPheno 3.2.1, Superiso 3.3 and SPheno by SARAH using the 
results normalized to the SM limit of each program. 

It is also possible to perform a check of self-consistency: the leading-order contribution has to be finite 
which leads to non-trivial relations among all wave and penguin diagrams and the involved vertices. There- 
fore, we can check these relations numerically by varying the renormalization scale used in all loop integrals. 
The dependence on this scale should cancel and the branching ratios should stay constant. This is shown in 
Figure [2] while single contributions can change by several orders the sum of all is numerically very stable. 
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Figure 1: The top figure: R = BR(B° /i+/i~)/BR(B° -> fJ, + (i~)sM for the constrained MSSM and a 
variation of m . The other parameters were set to Mxi% = 140 GeV, tan/3 = 10, /i > 0. In the middle m 
and Mi/2 were varied simultaneously, while tan/3 = 30 was fixed. In the bottom figure we show log(i?) as a 
function of tan /3, while m = M1/2 = 150 GeV were kept fix. In all figures A — and p > was used. The 
color code is as follows: superiso 3.3 (dotted black), SPheno 3.2.1 (dashed red) and SPheno by SARAH 
(solid blue). 
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Figure 2: The figure shows | Yl Fa\ penguin and | ^Iwavc as well as the sum of both | Fa\- Penguin and 
wave contributions have opposite signs that interchange between Q = 10 2 GcV to Q = 10 3 GeV. 

3.4- N on- super symmetric models 

We have focused our discussion so far on SUSY models. However, even if SARAH is optimized for the 
study of SUSY models it is also able to handle non-SUSY models to some extent. The main drawback at the 
moment for non-SUSY models is that the RGEs can not be calculated because the results of Refs. [HJ [32] 
which are used by SARAH are not valid in this case. However, all other calculations like the ones for the 
vertices, mass matrices and self-energies don't use SUSY properties and therefore apply to any model. 
Hence, it is also possible to generate SPheno code for these models which calculates B® d — > 11. The main 
difference in the calculation comes from the missing possibility to calculate the RGEs: the user has to provide 
numerical values for all parameters at the considered scale which then enter the calculation. We note that 
in order to fully support non-supersymmctric models with SARAH the calculation of the corresponding RGEs 
at 2-loop level will be included in SARAH in the future [43 . 

4. Conclusion 

We have presented a model independent implementation of the flavor violating decays B® d —> li in SARAH 
and SPheno. Our approach provides the possibility to generate source code which performs a full one-loop 
calculation of these observables for any model which can be implemented in SARAH. Therefore, it takes care 
of the necessity to confront many BSM models in the future with the increasing constraints coming from 
the measurements of 5° d —> 11 at the LHC. 
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Appendix A. Conventions 



Appendix A.l. Passarino-Veltman integrals 

We use in the following the conventions of [44] for the Passarino-Veltman integrals. All Wilson coefficients 
appearing in the following can be expressed by the integrals 

'Q 2 " 



B (0,x,y) 
A 

B 1 (x,y,z) 
C (x,y,z) 

Coo (*, 0,0,2, V, z) 

D (x,y,z,t) 
D 00 (x,y,z,t) 



A + f + In 
2 



V 



In 

x — y \x 



(!) 



4-D 



J E + fog 47T 



1/ ^B Q (x,y,z)-B (Q,y,z) 1 . . 
2^ z ~y) ~ ^B (x,y,z) 



1 

y- z 



y , y z , z 

log log - 

x — y x x — z x 



1 tz + xy — xz — yz + zJ 
4 + 2t 



C (0,0,t,x,z,y) 



+—[(t-x-y + 2z) B (t, x, y) + (x- z)B o (0, x, z) + (y - z)B o (0, y, z)] 



z log - 



t\og- r 



(y - x)(y - z)(y - 1) (z - x)(z - y)(z - t) (t - x)(t - y)(t - z) 
y 2 iog| z2l °9^ t 2 log |: 



(A.l) 
(A.2) 
(A.3) 
(A.4) 

(A.5) 

(A.6) 
(A.7) 



{y - x){y - z){y - t) (z - x)(z - y)(z - t) (t - x)(t - y)(t - z) _ 

Note, the conventions of ref. [H] (Pierce, Bagger [PB]) are different than those presented in ref. [3D] (Dedes, 
Rosiek, Tanedo [DRT]). The box integrals are related by 



D 
D Q0 



D 



D 



(PB) 


(PB) 
27 



-D, 



{DRT) 







D 



(DRT) 



(A.8a) 
(A.8b) 



Appendix A.2. Parametrization of vertices 

We are going to express the amplitude in the following in terms of generic vertices. For this purpose, we 
parametrize a vertex between two fermions and one vector or scalar respectively as 

Ga1u.Pl + Gb^Pr , (A.9) 

G a P l +GbPr. (A.10) 

Pl,r — |(1 T 7 5 ) ar e the polarization operators. In addition, for the vertex between three vector bosons 
and the one between one vector boson and two scalars the conventions are as follows 

Gvvv ■ (g^( k 2 - ki) P + 9u P {kz - k 2 ) f _ l + gp^h - k 3 )u) , (A. 11) 

Gssv ■ (h - fc 2 ) M . (A.12) 

Here, ki are the (ingoing) momenta of the external particles. 

Appendix B. Generic amplitudes 

We present in the following the expressions for the generic amplitudes obtained with FeynArts and 
FormCalc. All coefficients that are not explicitly listed are zero. Furthermore, the Wilson coefficients are 
left-right symmetric, i.e. 

Cxrr — Cxll{L <h> R) , Cxrl = Cxlr{L ^ R), (B-l) 

with X — S, V and where (L -H> R) means that the coefficients of the left and right polarization part of each 
vertex have to be interchanged. 
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Appendix B.l. Tree Level Contributions 




Figure B.3: Tree level diagrams with vertex numbering 



Since in models beyond the MSSM it might be possible that B° s — > il is already possible at tree-level. 



This is for instance the case for trilinear i?-parity violation [33J. The generic diagrams are given in Figure B.3 
The chiral vertices are parametrized as in eqs. ( A. 9 )-( A. 10 ) with A = 1, B = 2 for vertex 1 and A — 3, B = 4 



for vertex 2. Using these conventions, the corresponding contributions to the Wilson coefficients read 



r<( a ) 

U SLL 


1C 2 C1G3 


r (b) 
^VLL 


= 167T — = , 

My - S 


r (c) 


TP 2 -G1G3 

2(M| - t) 


r {d) 
SLR. 


= i6 *V-t« 


Me) 
^SLL 


= '"^(Ml-u) 


r (f) 
'-'SLR 


My — U 



Here, s, t and u are the usual Mandelstam variables. 
Appendix B.2. Wave Contributions 



WLR - 167T M , 



G1G4 

S 

G1G4 



Ml 



- s 

p" _|,,--' -G^Gs 
VLR ~ 2 (M| - t) 



C<* = 16tt , , , 
yii Ml 



G1G3 



C 



(e) 



(/) 

KM 



= 16tt 2 
16tt 2 ~ 



t 

G2G3 
2(M| - «) 
G1G4 



M 2 



(B.2) 
(B.3) 
(B.4) 
(B.5) 
(B.6) 
(B.7) 



The generic wave diagrams are given in Figure B.4 The internal quark which attaches to the vector or 
scalar propagator has generation index n. Couplings that depend on n carry it as an additional index. The 
chiral vertices are parametrized as in eqs. (A.9)-(A.10) with A = 1,B — 2 for vertex 1, A = 3, B = 4 for 




Figure B.4: Generic wave diagrams. For every diag: 
to the other external quark. 




there is a crossed version, where the loop attaches 
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Figure B.5: Generic wave diagram vertex numbering 



vertex 2, A — 5, B = 6 for vertex 3 and A = 7, B = 8 for vertex 4, see also Figure [BT5] for the numbering of 
the vertices. If a vertex is labelled 3' for instance, the corresponding couplings are G' 5 , G' 6 . Furthermore, we 
define the following abbreviations: 

fsi = a 1 2 (-MFidG^mn + G 2 G 4n m l )5< i 4) + (G 2 G 3n m n m t + G 1 G 4n m^)B^) , (B.8) 

fs2 = TO2 ^ m2 ( y M F (G 2n G i m j + G ln G 3 m n )B [ j) - (G 2n G 3 m 2 J + G ln G imj m n )B { ^ , (B.9) 

J n 

fs2 = 2 1 = (M F (Gi n G imj + G 2n G imn )B[p - (G ln G 4 m| + G 2n G 3 m 3 m n )B { / } ) , (B.10) 
M = — 2 (2M F (G 1 G 4 „m„ + G 2 G 3n mi)B^ + (G 2 G in m n mi + dG^m^B^) , (B.ll) 
fv2 = 2 l _ 2 (2M F (G 2 „G 3 m 3 + G ln G A rn n )B^ + (G 2 „G 4 m| + d^m^)^) , (B.12) 
M = m2 l _ m2 (2M F (G ln G im3 + G 2n G 3 m n )B { f) + {G ln G 3 m 2 } + G 2n G imj m n )B{^ . (B.13) 
The mi,rrij are the quark masses. Using these conventions, the contributions to the Wilson coefficients are 

C SLL = T72 1 ( G 5nfsi + G' 5n f S2 ) (B.14) 

M SQ - S 

c vll = .J? 7 (-G 5n fsi - G' 5n .fe) ->• M2 7 G 5 G x G 4 -Bi(0, M f , M s ) (B.15) 

m VQ ~ S V ; IVI VQ ~~ S 

C SLL = T72 (G 5n /vi - G' bn f V2 ) (B.16) 
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c vll ~ M 2 ~ : (-Gsnfvi + G'snfvi) -> A/f2 3— G 5 G 1 G 3 B 1 (0, M F , My) 

CfiLL = 2(Af 2 -t) ( G5nGr f S1 + G 'hn G ' 7 fs2) 

C VLR = 2(M| -t) ( G5 " G8 ^ S1 + G 6r l G 7/s2) 

CsLfl = » f2 + , ( G 5nG S fsi + G 6 „G' 7 / S2 ) 

C VLL = M 2 _ t ( G 5«G 7 /si + G' 5n G' 7 f S2 ^J 
r (f) 

^VLR — 772 7 I Lr 5n<^8.IVl - < ^6n Lr 7JV2 

r (h) 
SLR 

VLL ~ M 2 — t y^ 5 ^ 7 ^ 1 + U 5n Lr 7JV2 J 

c sll = 2(M 2 - u) ( G5 " G7 ^ S1 + G' 5n G' 7 .fs2) 
— (g^uGs/si + Gg„Gy/ss 



GyLL 



G { slr - 772 7 (G 6n G$f S i + G' 6n G' 8 fs2) 



VLR 
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Ko- 




-1 




M 2 S0 - 
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u 
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u 


-1 




M 2 S0 - 


u 


-4 




^0- 


u 


-2 




Ko- 


u 



C SLL - Ml -U V^W/Vl _ G 5n Cj 7/V2 ) 

C VLL = 772 7 (-G 5n G s f V i + Gg„G' 7 /y 2 ) 

m so u \ / 

SLR ~ M 2 - U \ L ' 6nLrS -' V1 ~ Lr 6n Cr 8/V2 J 

GyLH = M 2 " _ " (G &n G 7 fvi - G' 5n G' 8 fv2^ 

Appendix B.3. Penguin Contributions 
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"I 

Figure B.6: Vertex number conventions for a representative penguin diagram 
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(a) 



(b) 



(c) 



(d) 







(e) 



(f) 



00 







(i) 



0) 



(k) 



(1) 



Figure B.7: Generic penguin diagrams 



Diagrams with scalar propagators have Cvxy = and those with vector propagators have Csxy = 0. 
The vertex number conventions are given in fig. |B.6| and all possible diagrams are depicted in Figure |B.7| 



The chiral vertices are parametrized as in eqs. (A.9)-(A.10| with A = 1,B = 2 for vertex 1, A = 3,B = 4 



for vertex 2 and A = 7, B = 8 for vertex 4. Vertex 3 can be a chiral vertex, in this case A = 5, B = 6 is 
used. Otherwise, we will denote it with a index 5 and give as additional subscript the kind of vertex. The 
contributions to the Wilson coefficients from these diagrams read 



C 



(a) _ 
SLL " 



C 



r (b) 
VLL 



c 



(b) 

VLB, 



C 



C 



C 



(c) 
SLL 
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SLR 
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VLL 
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(a) _ 
SLR — 



M 2 S0 - 


s 
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M 2 S0 - 


— 
s 
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Mvo- 
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1 




M 2 - 

M S0 


— 
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M 2 - 

M S0 


— 
s 


2 




Ko 





G!G 3 G 7 ( G 6 B ( a b) + (G 5 M F1 M F2 + G 6 M|)C^ 



(a,b) 



y 



ddGs ( G 6 B ( a - b) + (-G 5 M F1 M F2 + G 6 M|)G^' 



GiG^G5 y sssGiM F C^' d ^ 
GiGsG5 t sssGsM F C c ' d 

104^5,55^^7^00 



(a,6) 



-,(a,6) 
y 00 



(B.34) 
(B.35) 
(B.36) 
(B.37) 
(B.38) 
(B.39) 
(B.40) 
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(0 

SLL 

w 

SLR 
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SLR 

eo 

(i) 
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Mso 



lLr4CT 5i S5y(_T 8 O 00 



-G1G4G7 (G 5 4 e ' /} + (G 6 M F1 M F2 + G 5 M V )C ( 



-G!G 4 G 8 [G 5 B ( eJ) + (G 6 M F1 M F2 + G 5 M v )C y 
-dGaGy (g^ J) + (-G 6 M F1 M F2 + G 5 M V )C^ J) - 2G 5 G, 

i \ 

-GiG 3 G 8 (g^ J) + (~G 6 M F1 M F2 + G 5 A/ 2 )G^ eJ) - 2G 5 G, 



Mvo 



Mvo 



00 

(e,f) 
00 



M| - S 



M| - 8 



GiG^G5^svvG'jM F C [ 



(g,h) 



GiGiG^^svvGg,M F C { 



(g,h) 




GiGzGayvvGi 



GiG^G^yvvGs 



(4' 
(4' 



(9,M 




M f Cq 9 ^ + 2C^ h) ) 



M F C { t h) + 2C { 9 h) ) 



M 



G 1 G 3 G^ssvG 7 [B^' + M F C' 



so 



j(»-0 , Mln^~ l ) 



— 2—— GiG 3 G 5 ,55yG 8 f_B l ^+M|-Gq 1 
M so s v 

-GiGiG$,ssvG-{ (^Bq ^ + M f Cq ^ 



M| - , 



M| - 



GiG4G5 j ssyG 8 ( f^ 1 ^ + M F C\ 



f 2 n {i-T) 



M V0 



-GiG4Gs t svv G 7 M F C K Q 
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-GiG4G5 i svvG 8 Mi?G 



^Vll — ~j^2 — ^~^C!iG3G 5 ^svvG 7 M F C 
1 ... ,-. „ 



(0 _ 

VLR 



M vo -s 



GiGzGz^svvGsMpCq 



Here, the arguments of the Passarino-Veltman integrals are as follows, with s — Mg : 



C ( x' b) =C x (s, 0,0, M F2 ,M F1 ,M t 



si 



C 



(c,d) 



= Cx(0, Sl 0,M F ,M 2 sll M 2 S2 ) 
C^ f) = C x (s,0,0,M F2 ,M F1 ,M v ) 
C { l h) = C x (0,s,0,M F ,M vl> M V2 ) 
c£~° = C x (0, s, 0, M F ,Mg, M v ) 



B^' b) =B x (s,M F1 ,M F2 ) 



B^ }) =B x {s,M 2 F1 ,M 2 F2 ) 



B 



{g,h) 



B x (s,M vl ,M 2 2 ) 



B x ~ l) =B x (s,M 2 ,M 2 ) 



(B.41) 
(B.42) 
(B.43) 
(B.44) 
(B.45) 
(B.46) 
(B.47) 
(B.48) 
(B.49) 
(B.50) 
(B.51) 
(B.52) 
(B.53) 
(B.54) 
(B.55) 
(B.56) 
(B.57) 



(B.58) 
(B.59) 
(B.60) 
(B.61) 
(B.62) 



Appendix B.4- Box Contributions 

The vertex number conventions for boxes are shown in figs. |B.8j while all possible generic diagrams are 
given in Figures B.9 and B.10 All vertices are chiral and they are parametrized as in eqs. (A.9)-(A.10) with 
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Figure B.8: Vertex number conventions for a set of representative box diagrams 



A = 1, B = 2 for vertex 1, A = 3, B = 4 for vertex 2, A = 5, £? = 6 for vertex 3 and A = 7, B = 8 for vertex 
4. If there arc two particles of equal type in a loop (say, two fermions), the one between vertices 1 and 2 (2 
or 3) will be labelled Fl and the other one will be F2. The contributions to the different Wilson coefficients 
read 

4ll = -G 1 G 3 G 5 G 7 M F1 M F2 ■ D^~ c) 

^slr = —GiGzGqG%MfiMf2 ■ Dy a c ' 

/~ii a ) s~i si s~i s~t n( a ~ c ) 

°VLL — -G 2 G 3 G 6 G7 ' ^00 

z^C a ) r* r< r* r* r>( a ~ c ) 

°VLiJ — -G 2 G 3 G 5 G 8 ' -^00 

4lx = -G 1 G^G b G 7 M F1 M F2 ■ ^ a ~ c) 

^slr = —GiGzGqG%M F iM F 2 ■ Dy a c ' 
/-<(fr) n n r 1 n( a— c ) 

n n c 1 r 1 n( a— c ) 

°VLfl — Cr 2 (_x 3 (_r 6 (jr 7 • -^00 

GgLL = -GiGsG^GrMpiM^D^ ^ 

C ( sl R = -2G 1 G 4 G 5 G 8 4 a - c) 

= -^G 2 G 4 G 5 G 7 M F1 M F2 4°- C) 

°VLi? - -G 2 G 3 G 5 G 8 -L» 00 
4ll = ^G^G^M, • D^ f) 

°SLi? - ^GiG 3 G 6 G 8 • -^oo 

r* r< r>(^ — /) 

VLL — — G 2 G 3 G 6 G7 • i^ 00 

Cvifl = ^G^G^M, • D { d ~ f) 
Csll = ^GiG 3 G5GyMi?iMF 2 • D ( d - f) 
^slr = 2GiG 3 G 6 G 8 • D ( Q q f) 
Gvll = ^2G2GzG^G%MfxMf2 ■ -Df) d ^ 

/^<( e ) s~t /~i fi fi r)(d—f) 

V LR ~ ( - T 2<-r3 ( - r 6 ( - r 7 ' u qo 

g'sll ~ 2GiG3GsGtMfiMf2D\) ^ 
c slr = -2GiG 4 G 5 G 8 L)^g }) 
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Figure B.9: Generic box diagrams I 



^vll — G 2 GiG^G 7 D m (B.85) 

C vIr = \G 2 G 3 G 5 G 8 M F1 M F2 D ( d - f) (B.86) 

41 = 2G 1 G 3 G e G 7 (Ct l) + M^D^ 21><g- i >) (B.87) 

cf LR = 2G 1 G 3 G 5 G & (C^ + M F1 D^ (B.88) 

C ( V 9) LL = G 2 G 3 G 5 G 7 M F1 M F2 D { a 9 - t) (B.89) 

GvIr = G 2 G 3 G 6 G 8 M F1 M F2 D { a 9 - l) (B.90) 

G ( sL = -4G 1 G 3 G 5 G 7 ^ 4) (B.91) 

4ifl = -4G 1 G 3 G 6 G 8 ^ 4) (B.92) 

C { y ] LL = G 2 G 3 G 5 G 8 M F1 M F2 D^- i} (B.93) 



1G 




(v) (w) (x) 

Figure B.10: Generic box diagrams II 



Cvlr = G 2 G 3 G 6 G 7 M F1 M F2 D { () g - t) (B.94) 

Cf LL = -G1G3G5G7 (C^ + MlD^ - 8D&-») (B.95) 

Cslfl = 2G 1 G 3 G 5 G 8 M F1 M F2 D { g - t) (B.96) 

C { Al = G 2 G 3 G 5 G 7 (C^ + MlDt i] 2D^) (B.97) 

Cvlr = G 2 G i G b G % M F1 M F2 D ( Q 9 - i) (B.98) 

C^ LL = 2G 2 G 3 G 5 G 7 (Ct l) + M 2 F1 Dt l) - 2^"°) (B.99) 

CsIr = 2G 2 G 3 G e G 8 (c^~ l) + M F1 D^~ 1 ^ - 2D^~ l) ) (B.100) 

d*\ L = GiG 3 G 6 G 7 M F1 M F2 D^~ 1 ^ (B.101) 

CvIr = GiG 3 Gc,G s M F1 M F2 D^~^ (B.102) 
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CsL = -4G 2 G 3 G 5 GsD ( >- l> (B.103) 

G { sL = -4G 2 G 3 G 6 G 7 Dg- l) (B.104) 

C { y ] LL = G 1 G 3 G 5 G 7 M F1 M F2 D i a J - l) (B.105) 

GvIr = G 1 G 3 G 6 G 8 M F1 M F2 D ( ^- l) (B.106) 

4L = -G 1 G 3 G 5 G 8 (ct l) + M 2 v D ( r l) - 8Dg~^) (B.107) 

Cf LR = 2G 1 G 4 G 5 G 7 M F1 M F2 D { j ~ l) (B.108) 

C$> LL = G 2 G 4 G 5 G 8 (C { ] - 1) + M v D^ l) - 2D&~ 1) ) (B.109) 

Cvlr = G 2 G 3 G 5 G 7 M F1 M F2 D^- [) (B.110) 

4ll = -GiG 3 G 6 G 7 ^ ro - o) + MlD { m - o) ^(UG^GeG, + 3G 2 G 4 G 5 G 8 )D^~ ^ (B.lll) 

4la = -2G 1 G 3 G 5 G 8 M i , 1 M F2 ^ m -° ) (B.112) 

GyL L = G 2 G 3 G 5 G 7 M F1 M F2 D ( a m - o) (B.113) 

= -G 2 G 3 G 6 G 8 (G< m - o) + M|i?( m - o) - 2^ - o) ) (B.114) 

4ll = 8GiG 3 G 6 G 8 L'o™ _0 ' 1 (B.115) 

4™Li? = 2G 1 G 3 G5G 7 M F1 M F2 L>' m " o) (B.116) 

G^ L - -2G 2 G 3 G 6 G 7 ^™- o) (B.117) 

Cvifl = G 2 G 3 G 5 G 8 M F1 M F2 ^ m - o) (B.118) 

G^ L = -^(13GiG 3 G 5 G 7 + 3G 2 G 4 G 6 G 8 )^™- o) (B.119) 

c sIr = -2G 1 G 4 G 5 G 8 M F1 M F2 ^ m_o) (B.120) 

Cy\ L = G 2 G 4 G 5 G 7 M F1 M F2 D ( ™- o) (B.121) 

CvL = 2G 2 G 3 G 5 G 8J D^- o) (B.122) 

4ll = -G 2 G 3 G 5 G 7 + M*D<T r) \(13G 2 G 3 G 5 G 7 + 3G 1 G 4 G 6 G 8 )£>^- r) ) (B-123) 

G^Lfl - -2G 2 G 3 G 6 G 8 M F1 M F2 D ( a p - r) (B.124) 

c££ L = GiG 3 G 6 G 7 M F1 M F2 ^ p_r) (B.125) 

4 P L = -GiG 3 G 5 G 8 (G( p - r) + M^t r » - 2D<£" r) ) (B-126) 

4ll = 8G!G 3 G 5 G 7J D^ r) (B.127) 

C^ LF = 2G 1 G 3 G e G 8 M F1 M F2 D ( J" r) (B.128) 

4ll = -2G 2 G 3 G 5 G 8 ^ r) (B.129) 

Gvlr = G 2 G 3 G 6 G 7 M F1 M F2 D^ p - r) (B.130) 

4ll = -^(13G 2 G 4 G 5 G 7 + 3G!G 3 G 6 G 8 )£)^- r) (B.131) 

C slr = —^G2GzG^G^,MfiMf2^{) + ^(G^G^G^GV — GiG?,G§G%)D m (B.132) 

= GiG 3 G 5 G 7 M F1 M F2 D^ p ~ r ^ (B.133) 

C^Lj = 2G 1 G 4 G 5 G 8J D^ r) (B.134) 
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C ( S S [ L = -4G 2 G 3 G e G 7 M F1 M F2 D ( s u> (B.135) 

4L = -4G 2 G 3 G 5 G 8 M F1 M F2 D { a s - u) (B.136) 

C ( vIl = -4GiG 3 G 5 G 7 (Ct u) + M F1 Dt u) 3£&- u) ) (B.137) 

4L = -^GsGeGs (g^ + M F1 D { Q s ~ u) ) (B.138) 

42l = -4G 2 G 3 G 5 G S M F1 M F2 D { S - U) (B.139) 

4L = -4G 2 G 3 G 6 G 7 M F1 M F2 D ( Q S - U) (B.140) 

C® LL = WGiGsGaGrD^ (B.141) 

4L = 4G 1 G 3 G 6 G 8J D^" ) (B.142) 

c££ L = -4G 2 G4G5G7M F1 M F2 L>^ u) (B.143) 

4i« = SG 2 G 3 G 5 G S D { S - U) (B.144) 

G^x - 16GiG 3 GbG7^'o" u) ( b - 145 ) 

C ( y ] LR = 2G 1 G 4 G 5 G 8 M F1 M F2 ^ s - l,) (B.146) 

c sll = $G 2 G 3 G 6 G 7 M F1 M F2 D ( J- X) (B.147) 

Cslr = 8G 2 G 3 G 5 G 8 (g^" x) + M vl D^ v ~ x) ) (B.148) 

4t L = -4G!G3G5G 7 (c^~ x) + M V1 D^ 3D&- X >) (B.149) 

C ( y ] LR = 2G 1 G 3 G 6 G & M F1 M F2 D^- X) (B.150) 

C^ L = ?,GxG 3 G % G % M fi M F2 D { q~ x) (B.151) 

4ifl = 32G 1 G 3 G 5 Grf - x) (B.152) 

G^l = -2G 2 G 3 G 6 G 7 M F1 M F2 D ( a v - x) (B.153) 

4^ - 4G 2 G 3 G 5 G 8 ^- a:) (B.154) 

C ( s x l L = -4GiG 4 G 5 G 8 M F1 M F2 £)^" x) (B.155) 

4la = -8G 1 G 3 G 5 G 7 (G^- a;) + M 2 2 D^ x) 3D^) (B.156) 

Gy\ L = -2G 2 G 3 G 5 G 8 M F1 M F2 ^- a:) (B.157) 

4L = -4G 2 G 4 G 5 G 7 (G^- a;) + M^Z^') (B.158) 

The arguments of the loop functions for the different amplitudes are 

4T C) = D X {M 2 F1 ,M 2 F2 ,M 2 S1 ,M 2 S2 ) (B.159) 

D%~ f) = D x (M 2 F11 M 2 F2l M 2 sll M 2 S2 ) (B.160) 

C { r i] = C x (0 3 ,M 2 F2 ,M 2 ,M 2 ) D^- l) = D X (M 2 F1 ,M 2 F2 ,M 2 ,M 2 S ) (B.161) 

4?"° - C X (6 3 ,M 2 2 ,M 2 ,M 2 ) D ( f l) = D x (M 2 F11 M 2 F2l M 2 ,M 2 ) (B.162) 

C , £ n ~° ) = C x (0 3 ,M 2 F2 ,M 2 F1 ,M 2 ) D x n - o) = D X (M 2 F2 ,M 2 F1 ,M 2 ,M V ) (B.163) 

4T r) =Gx(0 3 ,M| 2 ,M F1 ,M|) =D x (M 2 F2 ,M 2 F1 ,M Vl M 2 s ) (B.164) 

G^ u) = Gx(0 3 ,M F2 ,M2 1; M2 2 ) = D X (M 2 F1 ,M 2 2 ,M V1 ,M 2 2 ) (B.165) 

c£~ x) = C x (d 3 ,M 2 2 ,M F1 ,M 2 2 ) D%- x) = D X (M F2 ,M 2 F1 ,M V1 ,M 2 2 ) (B.166) 
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